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ABSTRACT 

We report the discovery of 13 confirmed two-image quasar lenses from a systematic search 
for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey 
(BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search 
(SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular 
separations (< 2") confirmed from our spectroscopy, which are either projected pairs, physi¬ 
cal binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. 
The newly discovered quasar lens system, SDSS J1452-1-4224 at 4.8 is one of the high¬ 
est redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens 
candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the 
lens sample presented here does not have a well-defined selection function. 

Key words: gravitational lensing: strong - quasars: general methods: observational - meth¬ 
ods: statistical 



2 More et al. 


1 INTRODUCTION 


Gravitationally lensed quasars provide a unique tool to study the 
Universe. In particular, the time-variable nature of quasars enable 
us to measure the arrival time difference between quasar multiple 
images, which may be a powerful probe of cosmology (e.g., |Refs^ 
|dal||1964{ ISchechter et al.|ri997[ |Treu et 5I]|2013| l. Applications 
of gravitationally lensed quasars such as time delay cosmography 
(e.g.,|Coles|2008[[Suyu et al.|2010^, constraints on the quasar lumi 


nosity function (e.g.,|Comerford et al.|2002 Richards et al. 


20061, 


20121 , 


constraints on dark energy (e.g., |Kochanek|199^ Oguri et al. 
and study of the host galaxy and/or black hole properties (e.g., |Peng| 
|et al.|2006{|Mosquera et al.|2013[[Rusu et al.|2015| l are still limited 
by small number statistics, indicating the importance of finding new 
gravitationally lensed quasars in wide-held surveys. 

The SDSS Quasar Lens Search (SQLS; |Oguri et al.||2006| 
|2008[[2OT^|Inada et al.|2008||2010l|2012[ l considerably advanced 
the held by discovering nearly 50 new quasar lenses. The SQLS 
was based upon the large sample of spectroscopically conhrmed 
quasars in the Sloan Digital Sky Survey - I/II (SDSS; [York et al.| 
|2000| >. Candidates were hrst selected from the ~ 100,000 spectro¬ 
scopic quasars and their SDSS images were then examined to iden¬ 
tify quasars with extended morphology or with nearby companion 
objects of similar colour. The lens candidates were then conhrmed 
with various facilities to construct a secure sample of gravitation¬ 
ally lensed quasars. Thanks to the well-studied selection function 
( [Oguri et al.|2006| >, a sub-sample of the SQLS quasar lenses was 
used to place statistical limits on dark energy and the evolution of 
massive galaxies that act as lenses ( [Oguri et al.|2012^ . 

The Baryon Oscillation Spectroscopic Survey (BOSS; [Daw-[ 
[son et al.|2013[ l of the SDSS-III jEisenstein et al.|2011[ l is a new 
spectroscopic survey to measure the expansion rate of the Universe. 
In addition to an extensive spectroscopic survey of 2 ~ 0.5 galax¬ 
ies, BOSS has obtained spectra of ~ 300,000 quasars ( [Alam et^ 
[2015[ l with over 180000 in the redshift range 2.15 < 2 < 4 in order 
to detect baryon acoustic oscillation signatures in the Lya absorp¬ 
tion features of their spectra. The large number of BOSS quasars 
suggests that it should contain a large number of gravitationally 
lensed quasars. However, the lensing rate of the BOSS quasar sam¬ 
ple is likely smaller than that of the SDSS-II quasar sample be¬ 
cause BOSS largely targets point-like quasar candidates in order to 
increase the survey efficiency. In this case, lensed quasars where 
the images are not individually resolved would appear extended 
and would tend to be excluded by the BOSS selection function. 
On the other hand, the typical BOSS quasar is at a higher redshift 
(2 ~ 2 — 3 [Ross et al.[[2012l [Paris et al.[[2014[ l than the typical 
SDSS-I/II quasars (2 ~ 1 — 2) and this higher average redshift may 
partly compensate for the lower lensing rate due to the morpholog¬ 
ical target selection. 

We thus started a new survey for gravitationally lensed quasars 
from BOSS called the BOSS Quasar Lens Survey (BQLS). The 
basic strategy follows that of the SQLS in that we select quasar 
lens candidates by examining the SDSS images of BOSS quasars 
and conduct additional observations of any promising quasar lens 
candidates. In this paper, we present 13 confirmed quasar lenses 
from the BQLS. In addition, we also identify 11 quasar close pairs 
with not quite identical but similar spectra. Some of these may be 
true gravitational lenses but could not be confirmed due to lack of 
detection of a lens galaxy or spectra being slightly different due to 
effects such as microlensing and dust extinction. 

This paper is organized as follows. In Sectionj^ we describe 
the selection of our lens candidates. In Section we describe the 


additional observations undertaken for a detailed examination of 
the most promising lens candidates. The mass modelling of the 
lenses is described in Section [^ We compare the confirmed lenses 
to the SQLS lens sample and present our general findings in Sec- 
tionj^ Sectionj^presents our summary and conclusions. 


2 CANDIDATE SELECTION 

We selected candidate quasar lens systems from the imaging and 
spectroscopic data of SDSS-III/BOSS. The SDSS uses a dedicated 
2.5-meter telescope at Apache Point Observatory in New Mex¬ 
ico ( [Gunn et al.|2006^ . The SDSS imaging data uses five (ugriz) 
broadband filters and has a limiting magnitude of r ~ 22.5, which 
were used to select various targets including quasars for spectro¬ 
scopic observations jRoss et aL][2012^ . The SDSS-III/BOSS sur¬ 
vey uses an upgraded multi-object fiber spectrograph which can 
simultaneously take spectra of 1000 objects with a resolving power 
R ~ 2000 and a wavelength coverage of 3600A < A < lOdOOA 
( [Smee et al.[2013j l. 

We largely followed the procedures established in the SQLS 
for the lens candidate selection process. In the SQLS, however, 
there were also image separation and flux ratio limits for construct¬ 
ing a complete lens sample for statistical studies. In SDSS-III, con¬ 
structing a well-defined statistical sample is much more challeng¬ 
ing because of the complex nature of the quasar target selection 
( [Ross et al.|2012^ . Thus, we first adopted a simple approach of ap¬ 
plying the SQLS-like selection criteria and then choosing the best 
lens candidates by visually inspecting the SDSS images of all of 
the resulting lens candidates. Below we describe these procedures 
in more detail. 

The selection method for the underlying SDSS-III/BOSS 
spectroscopic quasar sample is presented in [Ross et ^ \2012\ 
and [Paris et al.[ j2014^ . Currently all of the BOSS data 
are public, and the quasar catalogue we used is available 
at http://www.sdss.org/drl2/algorithms/boss-drl2-quasar-catalog/. 
From the spectroscopic quasar catalog, we selected quasar lens can¬ 
didates using a method similar to SQLS, in which morphological 
and colour criteria were used to select lens candidates from the ob¬ 
ject catalogue. The full details of the SQLS selection algorithm are 
presented in [Oguri et al.[ ( |2006) and [Inada et al.[p008^ . While the 
SQLS selection criteria were optimized for selecting low-redshift 
(2 < 2.2) quasar lens candidates, an extension of the selection al¬ 
gorithm to higher redshifts is discussed in [Inada et^Tj ( [2009[ l. In 
this paper, we adopted a modified version of these morphological 
and colour selection criteria as described below. 

The morphological selection approach identified spectro¬ 
scopic quasars that are poorly fit by the imaging Point Spread Func¬ 
tion (PSF). Fven though BOSS targeted point-like quasars, there 
are extended objects with a spectroscopic classification of a quasar 
in BOSS which were originally identified as galaxies and targeted 
as part of the galaxy follow-up programs or ancillary programs tar¬ 
geting AGNs. Therefore, we applied the morphological selection 
to find small image separation lensed quasars which appear as ex¬ 
tended quasar-like sources. Specifically, we used the following ob¬ 
ject parameters; type, which describes the classification of objects 
as stars (type = 6) or galaxies (type = 3); InlStar, which is 
the logarithm of the likelihood of objects being fitted by the PSF 
of the field; InlDeV, which is the logarithm of the likelihood of 
objects being fitted by the de Vaucouleurs profile. The parameters 
InlStar and InlDeV are available for each broad band. For 
quasars at 0.7 < 2 < 3 we selected lens candidates whose gri- 
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band object parameters satisfy InlStar < —1 for type = 3 and 
InlDeV < —15 and InlStar — InlDeV < —5 for type = 6. For 
quasars at z > 3, we used the same criteria as above but only for 
the r and i-bands. There were ~ 3000 candidates selected hy these 
criteria. 

For a given spectroscopic quasar, the colour selection method 
searches for nearby objects with similar colours out to a maximum 
image separation of A6 = 30”. For quasars at 0.7 < 2 < 3.4 we 
used all 5 bands and select objects with colour differences smaller 
than 0.1 for g — r and r — i and smaller than 0.2 for u — g and 
i — z, where the colours were computed using PSF magnitudes. 
For quasars at 2 > 3.4, we used only the griz-bat\d& and similar 
thresholds for the colour differences of 0.1 for r — i and 0.2 for 
g — r and i — z. We also required that the i-band magnitude dif¬ 
ference of the two objects was smaller than 1. When the angular 
separation of the two objects were small (< 375) we used relaxed 
colour criteria as the SDSS photometry might not be reliable (see 
|Oguri et al.|2006[ l. This colour selection method led to a sample of 
~ 2000 candidates. 

We then visually inspected all the morphological and colour 
selected lens candidates. The visual inspection was performed by 
one of the authors (MO) to define “good” lens candidates that have 
relatively high chance of being true lens systems. Here the good¬ 
ness of the candidates was judged by various factors including 
the presence of multiple components with quasar-like colours or 
nearby galaxies that might act as lenses. Most of the false positives 
removed in the visual inspection are either single quasars with con¬ 
tamination from bright nearby objects and obvious quasar-galaxy 
superpositions. Large-separation (> 5”) quasar lenses are much 
less common, and hence large-separation lens candidates are not 
regarded as good unless putative lens galaxies or clusters are seen 
in between the candidate quasar images. We finally selected ~ 80 
good candidates and an additional ~ 70 possible lens candidates 
from the morphological selection, and ~ 75 good candidates and 
~ 65 possible lens candidates from the colour selection. There is 
an overlap of ~20 good candidates between the morphology and 
colour selected samples. The number of lens candidates, requiring 
observational follow-up studies was reduced by a factor of ~ 20 
through the visual inspection. This suggests that our approach is 
an efficient approach to find the most promising lens candidates. 
Nevertheless, it is possible that we missed some true gravitational 
lens systems in the course of the visual inspection. For example, 
we are likely to miss quasars that are fainter and do not visually 
look star-like which includes some small image-separation lenses 
and possibly quadruply imaged systems. We may also have missed 
systems with atypical image configurations, odd colours or highly 
anomalous flux ratios. 

We note that |Dahle et'H^ j2013| |2015^ discovered a gravi¬ 
tationally lensed sextuple quasar with a maximum image separa¬ 
tion of 15'.'1, SDSS J2222-H2745, from the SDSS-III data. They 
searched for strong lens candidates by visually inspecting photo¬ 
metrically identified clusters of galaxies. Even though this lens sys¬ 
tem meets our image separation criteria, none of the lensed quasars 
from SDSS J2222-I-2745 were selected for spectroscopy by BOSS 
so this source is not in our parent sample. 

Some gravitationally lensed quasars have been discovered by 
identifying quasar spectral features superposed on the spectra of 
lens galaxies (e.g., [Johnston et al.|2003t [Bolton et al.| |2005[ |Mc-| 
[Greer et al.|2010^ although this technique of identifying compound 
objects spectroscopically has primarily discovered many examples 
of lensed star-forming galaxies (e.g., [Bolton et al.[[2006b^ . As a 
complementary approach to the SQLS selection technique, we also 


searched for lensed quasars using this technique. We identified in¬ 
stances of quasar-galaxy superpositions in the BOSS spectra with 
quasars at 2 < 5.5 and galaxies with redshifts lower than the 
quasars (J. Zinn et al. in prep.). The initial sample had 17 good 
candidates out of which 3 systems were followed-up spectroscop¬ 
ically as part of the BQLS sample. Currently, this search has pro¬ 
duced ~100 candidates with a “good” visual quality flag and will 
be presented in Zinn et al. (in prep.). 

The BOSS DR12 catalogue has some known lenses, several of 
which are in the SQLS. This known sample was excluded prior to 
the visual inspection of the candidates. In total, we identified about 
250 new lens candidates for further study. In the next section, we 
present initial follow-up results of the most promising quasar lens 
candidates. 


3 DATA: OBSERVATION, REDUCTION AND ANALYSIS 

Our follow-up observations consist of imaging and spectroscopic 
observations. We first attempt to obtain follow-up images for the 
good lens candidates, and conduct follow-up spectroscopic obser¬ 
vations only when the follow-up images show evidence for gravita¬ 
tional lensing. However, it was not possible to follow this strategy 
strictly due to observing constraints. In some cases we obtained 
follow-up spectra before we conducted imaging follow-up obser¬ 
vations. Among these systems, those with clear differences in the 
spectra were thus not observed with imaging. Since the spectro¬ 
scopic observations are very informative for confirming or ruling 
out the lensing hypotheses of the candidates, below we first present 
the complete spectroscopic follow-up results for the best 55 can¬ 
didates, and then present imaging follow-up results for confirmed 
gravitational lens systems. The full imaging follow-up results will 
be presented in an upcoming paper. 

3.1 Spectroscopy 

Table [T] summarizes our spectroscopic observations. We obtained 
spectra of five candidates with the Low Resolution Imaging Spec¬ 
trometer (LRIS; [Oke et al.|1995[ > on the Keck telescope on 2013 
September 7. We used the I'.'O slit, the 400/3400 grism for the blue 
channel and the 400/8500 grating for the red channel. The spectral 
resolution was R ~ 600 in the blue channel and R ~ 1300 in the 
red channel. The red channel was binned to give 0727 pixel“^. The 
unbinned blue channel has 07135 pixel” We also obtained spec¬ 
tra of 50 candidates with the Faint Object Camera and Spectrograph 
(FOCAS;[Kashikawa et al.|2002f on the Subaru telescope on 2014 
May 2, 2015 Feb 19, and 2015 June 19. In all of the FOCAS obser¬ 
vations we used the l70 slit, but adopted two different configura¬ 
tions depending on the quasar redshifts; one is the grism and filter 
set of 300B/L600 to cover the wavelength range of 3700 — 6000 A, 
and the other is 300B/SY47 4700 — 9100 A. Both datasets have a 
spectral resolution of i? ~ 400 — 500. The data were binned 2x2 
on the detector, providing a spatial resolution of 0721 pixel”^. In 
all of the observations, the long-slit was aligned to observe both pu¬ 
tative quasar images simultaneously. The data were reduced using 
standard IRAlQtasks. 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Table 1. Summary of spectroscopic observations. 


Object ID 

RA (J2000) 

Dec (J2000) 

Selection 

ZQSO 

Observation 

Exp. [s] 

Result 


SDSS J0033+2015 

00:33:37.59 

+20:15:38.2 

c 

2.701 

K, 2013 Sep (400/3400+400/8500) 

600 

quasar+star 


SDSS J0035+2659 

00:35:31.98 

+26:59:59.9 

M 

2.294 

S, 2015 Jun (300B/L600) 

1500 

different SED 


SDSS JOl 14+0722 

01:14:38.38 

+07:22:28.5 

MC 

1.828 

K, 2013 Sep (400/3400+400/8500) 

600 

quasar lens {zi ■ 

= 0.408) 

SDSS J0139+1908 

01:39:40.69 

+19:08:40.7 

MC 

3.095 

K, 2013 Sep (400/3400+400/8500) 

600 

different SED 


SDSS J0206+0440 

02:06:49.50 

+04:40:19.0 

C 

2.396 

K, 2013 Sep (400/3400+400/8500) 

600 

different SED 


SDSS J0213-0421 

02:13:22.86 

-04:21:34.3 

C 

1.910 

S, 2015 Feb (300B/L600) 

900 

quasar pair z = 

1.911&0.992 

SDSS J0256+0153 

02:56:40.76 

+01:53:29.3 

s 

2.600 

K, 2013 Sep (400/3400+400/8500) 

600 

quasar lens {zi ■ 

= 0.603) 

SDSS J0737+4825 

07:37:08.67 

+48:25:51.1 

M 

2.892 

S, 2015 Feb (300B/L600) 

1080 

quasar lens 


SDSS J0757+2150 

07:57:20.54 

+21:50:07.7 

MC 

2.128 

S, 2014 May (300B/L600) 

720 

quasar+galaxy at 2 : = 0.118 

SDSS J0818+0601 

08:18:30.46 

+06:01:38.0 

M 

2.381 

S, 2015 Feb (300B/L600) 

900 

quasar pair 2 : = 

2.359&2.361 

SDSS J0821+0735 

08:21:43.36 

+07:35:45.9 

C 

2.378 

S, 2015 Feb (300B/L600) 

1080 

quasar pair 2 : = 

2.383&2.387 

SDSS J0821+4542 

08:21:58.66 

+45:42:44.4 

M 

2.066 

S, 2014 May (300B/L600) 

900 

quasar lens {zi ■ 

= 0.349) 

SDSS J0826+4248 

08:26:52.99 

+42:48:17.8 

MS 

1.942 

S, 2014 May (300B/L600) 

900 

inconclusive 


SDSS J0847+1504 

08:47:21.96 

+15:04:50.2 

M 

3.280 

S, 2015 Feb (300B/SY47) 

900 

different SED 


SDSS J0921+2854 

09:21:15.47 

+28:54:44.3 

MC 

1.410 

S, 2015 Feb (300B/L600) 

720 

quasar lens {zi ■ 

= 0.445“) 

SDSS J0928+4332 

09:28:39.19 

+43:32:42.4 

MC 

3.698 

S, 2014 May (300B/SY47) 

900 

quasar pair 2 = 

3.694&2.995 

SDSS J0930+4614 

09:30:21.16 

+46:14:22.8 

M 

2.397 

S, 2015 Feb (300B/L600) 

900 

quasar pair 2 = 

2.393&2.394 

SDSS J0958+0744 

09:58:42.24 

+07:44:23.2 

M 

2.781 

S, 2014 May (300B/L600) 

900 

quasar+galaxy at 2 = 0.158 

SDSS JlOOl+0156 

10:01:40.95 

+01:56:43.1 

M 

2.202 

S, 2015 Feb (300B/L600) 

1080 

different SED 


SDSS J1043+4320 

10:43:24.87 

+43:20:49.4 

M 

2.234 

S, 2015 Jun (300B/L600) 

900 

quasar pair 2 = 

2.245&2.225 

SDSS J1108+4726 

11:08:28.25 

+47:26:21.6 

C 

3.627 

S, 2014 May (300B/SY47) 

900 

quasar+galaxy at 2 = 0.385 

SDSS J1124+5710 

11:24:55.24 

+57:10:56.5 

C 

2.312 

S, 2015 Feb (300B/L600) 

720 

quasar pair 2 = 

2.309&2.315 

SDSS J1125+5020 

11:25:42.60 

+50:20:35.5 

c 

1.066 

S, 2015 Jun (300B/L600) 

900 

quasar+galaxy at 2 = 0.385 

SDSS J1138+5254 

11:38:38.54 

+52:54:18.2 

M 

2.758 

S, 2015 Jun (300B/L600) 

1200 

different SED 


SDSS J1139-0014 

11:39:28.49 

-00:14:18.1 

c 

3.084 

S, 2014 May (300B/L600) 

720 

quasar+star 


SDSS J1152+2235 

11:52:10.58 

+22:35:19.2 

c 

2.902 

S, 2014 May (300B/L600) 

720 

quasar+star 


SDSS J1248+6104 

12:48:59.85 

+61:04:30.3 

c 

2.591 

S, 2015 Feb (300B/L600) 

1200 

quasar+star 


SDSS J1254+1857 

12:54:40.37 

+18:57:12.0 

MC 

1.717 

S, 2015 Feb (300B/L600) 

1200 

quasar lens {zi ■ 

= 0.555“) 

SDSS J1309+5617 

13:09:27.55 

+56:17:38.9 

C 

2.505 

S, 2015 Jun (300B/L600) 

1200 

quasar pair 2 = 

2.513&2.515 

SDSS J1319+5023 

13:19:26.06 

+50:23:05.5 

C 

2.301 

S, 2015 Jun (300B/L600) 

1080 

quasar+galaxy at 2 = 0.088 

SDSS J1322+3038 

13:22:04.81 

+30:38:38.2 

M 

2.233 

S, 2015 Jun (300B/L600) 

1200 

inconclusive 


SDSS J1330+3800 

13:30:07.34 

+38:00:42.4 

C 

2.254 

S, 2015 Feb (300B/L600) 

720 

quasar lens 


SDSS J1405+1350 

14:05:56.92 

+13:50:38.3 

C 

2.342 

S, 2015 Jun (300B/L600) 

1200 

quasar pair 2 = 

2.345&2.361 

SDSS J1412+5204 

14:12:10.18 

+52:04:23.3 

M 

2.952 

S, 2015 Feb (300B/L600) 

1200 

quasar lens 


SDSS J1429+2523 

14:29:38.20 

+25:23:43.4 

M 

3.904 

S, 2015 Feb (300B/SY47) 

900 

different SED 


SDSS J1442+4055 

14:42:54.79 

+40:55:35.6 

C 

2.575 

S, 2014 May (300B/L600) 

480 

quasar lens 


SDSS J1452+4224 

14:52:11.50 

+42:24:29.6 

M 

4.819 

S, 2015 Feb (300B/SY47) 

900 

quasar lens {zi ■ 

= 0.382) 

SDSS J1458-0202 

14:58:47.59 

-02:02:05.9 

M 

1.724 

S, 2014 May (300B/L600) 

1200 

quasar lens 


SDSS J1508+3037 

15:08:22.32 

+30:37:47.2 

M 

2.449 

S, 2015 Feb (300B/L600) 

1200 

inconclusive 


SDSS J1537+3014 

15:37:34.46 

+30:14:53.7 

MCS 

1.553 

S, 2014 May (300B/L600) 

900 

quasar lens ( 2 ^ = 

= 0.490“) 

SDSS J1548+2830 

15:48:50.76 

+28:30:14.3 

C 

3.227 

S, 2015 Feb (300B/SY47) 

900 

quasar pair 2 = 

3.208&1.487 

SDSS J1551+3303 

15:51:37.02 

+33:03:19.6 

c 

3.806 

S, 2015 Feb (300B/SY47) 

900 

different SED 


SDSS J1552+2401 

15:52:21.40 

+24:01:23.1 

MC 

3.674 

S, 2015 Feb (300B/SY47) 

720 

quasar+star 


SDSS J1554+2616 

15:54:11.53 

+26:16:35.7 

C 

2.321 

S, 2015 Feb (300B/L600) 

720 

quasar+star 


SDSS J1556+1731 

15:56:23.81 

+17:31:21.4 

C 

2.814 

S, 2014 May (300B/L600) 

720 

quasar+galaxy at 2 = 0.109 

SDSS J1600+3148 

16:00:28.78 

+31:48:31.7 

C 

2.331 

S, 2015 Jun (300B/L600) 

1080 

quasar+star 


SDSS J1611+0844 

16:11:05.64 

+08:44:35.4 

c 

4.548 

S, 2015 Feb (300B/SY47) 

900 

different SED 


SDSS J1627+5553 

16:27:16.69 

+55:53:37.0 

c 

4.072 

S, 2015 Feb (300B/SY47) 

720 

different SED 


SDSS J1712+2516 

17:12:32.14 

+25:16:24.6 

M 

2.604 

S, 2015 Feb (300B/L600) 

1080 

quasar+star 


SDSS J1715+2831 

17:15:17.32 

+28:31:29.5 

c 

2.023 

S, 2015 Feb (300B/L600) 

1080 

different SED 


SDSS J2146-0047 

21:46:46.03 

-00:47:44.3 

M 

2.381 

S, 2014 May (300B/L600) 

900 

quasar lens at 2 

= 0.799 

SDSS J2220-0050 

22:20:41.31 

-00:50:15.5 

MC 

2.602 

S, 2015 Jun (300B/L600) 

1080 

quasar+star 


SDSS J2238+1820 

22:38:51.80 

+18:20:38.5 

MC 

2.074 

S, 2015 Jun (300B/L600) 

720 

quasar+galaxy at 2 = 0.199 

SDSS J2245+2548 

22:45:55.75 

+25:48:35.4 

C 

2.992 

S, 2015 Jun (300B/L600) 

900 

quasar pair 2 = 

2.995&2.178 

SDSS J2314-0108 

23:14:24.50 

-01:08:58.7 

C 

3.455 

S, 2015 Feb (300B/SY47) 

1080 

different SED 



The ‘Selection’ column indicates the method(s) with ‘M’ for morphological selection, ‘C’ for colour selection, and ‘S’ for spectroscopic selection (see text 
for details). Several candidates were selected by multiple methods. ‘Observation’ indicates the telescope (S—Subaru or K—Keck), the date and setup of the 
spectroscopic observations (see also Sec. HD ‘Result’ indicates the conclusion from the follow-up spectroscopy as well as imaging observations described in 
Sec. |3.2| The 2 :qso are taken from the BOSS DR12 catalogue which has redshifts corrected after visual inspection of the spectra. The last column has 
redshifts from the follow-up spectra except when marked with “ which are taken from the BOSS spectrum. 
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Table 2. Summary of imaging observations of the confirmed lenses. 


Table 4. Sersic parameters of the lens galaxies. 


Object 

ID 

Instrument 

Observing Date 

Exposure 

(sec) 

Object 

ID 

Te 

n 

n 

e 

de 

(deg) 

J0114-t0722 

Tek2k 

2013 Dec 24 

1440 (V,I) 

J0114-t0722 

0.84±0.04 

2.13±0.24 

0.23±0.02 

73.3±6.2 

J0256-t0153 

Tek2k 

2013 Dec 24 

1440 (V,I) 

J0256-t0153 

1.16±0.06 

1.35±0.21 

0.33±0.03 

72.5±4.3 

J0737-t4825 

Tek2k 

2015 Feb 20 

1440 (V,I) 

J0737-t4825 

1.08±0.35 

[4.00] 

0.14±0.12 

-70.0±32.6 

J0821-t4542 

Tek2k 

2014 Feb 21 

2160 (V) 1980 (R) 

J0821-t4542 

0.14±0.02 

[4.00] 

[0.00] 

[0.0] 

J0921-t2854 

Tek2k 

2015 Feb 20 

1440 (V) 1080 (I) 

J0921-t2854 

0.70±0.01 

0.98±0.06 

0.10±0.02 

36.7±9.6 

J1254-tl857 

Tek2k 

2014 Feb 21 

3060 (V) 2880 (R) 

J1254-tl857 

0.73±0.14 

[4.00] 

0.38±0.12 

-16.1±14.5 

J1330-t3800 

Tek2k 

2014 Feb 21 

1440 (V) 2160 (R) 

113304-3800 

0.87±0.34 

[4.00] 

0.75±0.16 

-16.3±9.7 

J1412-t5204 

Tek2k 

2015 Feb 20 

1440 (V,I) 

114124-5204 

0.61±0.15 

[4.00] 

0.58±0.08 

15.5±7.4 

J1442-t4055 

OSMOS 

2013 Apr 08 

1200 (g) 2100 (r) 900 (i,z) 

J14424-4055 

0.99±0.10 

4.19±0.70 

0.21±0.03 

68.2±6.1 

J1452-t4224 

FOCAS 

2015 Jun 19 

270 (R) 360 (I) 

114524-4224 

1.09±0.01 

[4.00] 

0.30±0.01 

-23.8±1.1 

J1458-0202 

SOI 

2013 Feb 12 

420 (g,r,i) 

11458-0202 

2.13±0.24 

3.81±0.30 

0.36±0.02 

30.6±1.6 

J1537-t3014 

Tek2k 

2014 Feb 21 

1440 (V) 1800 (R) 

115374-3014 

3.57±2.09 

5.90±1.92 

0.37±0.04 

87.2±5.1 

J2146-0047 

Tek2k 

2012 Sep 17 

960 (V,I) 

12146-0047 

0.38±0.28 

[4.00] 

[0.00] 

[0.0] 


Tek2k is on the UH88 telescope, OSMOS is on the Hiltner telescope, 

FOCAS is on the Subaru telescope and SOI is at the SOAR telescope. ^he parameters n, e, and 6, denote the effective radius, Sersic index, 

ellipticity, and the position angle (measured counter-clockwise from North 
where North is up at 0 deg), respectively. See Table|^for the magnitudes. 


We found that thirteen out of the 55 quasar lens candidates 
contain two quasar components with identical redshifts and simi¬ 
larities in the Spectral Energy Distributions (SED). Together with 
the analysis of the imaging observations presented below, we con¬ 
clude that these 13 systems are gravitationally lensed quasars. Our 
spectroscopy also identified an additional 11 quasar pairs with ei¬ 
ther almost identical or slightly different redshifts. Properties of 
these 11 quasar pairs will be discussed in Sec. |3.3| and |5.3| 

Figure[T]presents the spectra for the 24 confirmed quasar pairs, 
including the 13 confirmed gravitationally lensed quasars. Eor five 
of the quasar lens systems, the lens galaxies are bright enough to 
be clearly detected in one of the quasar spectra. We successfully 
measured the lens redshift for these lens systems (see Table[TJ. For 
an additional three quasar systems, the lens galaxies are bright and 
their redshifts were derived from the SDSS-III spectra. The subse¬ 
quent spectra do not cover the red wavelength range where these 
lens galaxies can be detected easily. 

3.2 Imaging 

We imaged a sub-sample of the most promising candidates in 2 — 3 
bands with the goal of detecting the lens galaxy. At minimum, 
these imaging data should better resolve the small separation lensed 
quasar images to provide the basic constraints required for mass 
modelling. The majority of the lenses from the BQLS sample were 
observed with the Tektronix 2048 x 2048 CCD camera (Tek2k, with 
a pixel resolution of 0”22) on the UH88 telescope. A few of the 
candidates were observed with other telescopes. A summary of 
these observations is given in Table which includes the exposure 
times for each filter. We used the FOCAS instrument with 2x2 
binning on Subaru giving a pixel resolution of 0”208, the SOI in¬ 
strument (pixel resolution of O'.'IS) on the SOuthern Astrophysical 
Research (SOAR) telescope and the OSMOS instrument (pixel res¬ 
olution of 0”273) on the Hiltner telescope of the MDM observatory. 
We processed the data in IRAF by following standard procedures. 
The seeing of our imaging data is better than 1". 

In most cases, the imaging reveals two point-like components 
with a hint of extended emission close to one of those components 
indicating presence of a lens galaxy. We use GALFIT ( |Pen g et aT] 
|2002| [^10|> to model each of these components and to measure 


their fluxes and positions. The quasars are assumed to be point 
sources convolved with the PSF and the galaxy is modelled with a 
Sersic profile. The use of point sources for the quasar components 
is reasonable because all the quasars are at high redshifts 2 > 1. 
We use the nearby stars as a model for the PSF. We first fit models 
to the two quasars without the lens. Next, we fit models including 
the lens. Both modelling results are compared and the latter model 
is accepted if the and the residuals have improved. In the cases 
when it was difficult to fit the Sersic index of the galaxy’s light 
profile, we held it fixed at n = 4, since the de Vaucouleurs profile 
( |de Vaucouleurs|1948| ) is a good model for early-type galaxies. The 
modelling results for all of the 13 confirmed lenses are presented in 
Figure]^ For each lens, we show the reddest available image (left 
panel) where the contribution from the lens galaxy is easier to see 
than in the other bands, the GALFIT model which includes the lens 
galaxy (second panel from left) and the residuals from this model 
(third panel from left). For comparison, we also show the residuals 
from the model without the lens galaxy in the right-most panel. 

We performed relative photometry of the lenses by comparing 
them to nearby stars with known SDSS magnitudes. First, we 
used the software from Astrometry.net ( |Lang et al.||2010l > to 
determine the orientations of the images and the positions of 
the stars, except for SDSS J1458-0202 and SDSS J1452-I-4224 
where the astrometry was determined manually. For every lens 
in every band available from the imaging, we measured the 
fluxes of about a dozen stars using Sextractor ^Bertin &| 
|Amouts||1996^ . The zero points were determined by comparing 
the SDSS magnitudes to the fluxes of the stars. Owing to the 
differences between the SDSS filters and those used for the follow 
up observations, we had to apply appropriate filter conversions 
to the fluxes. We adopted the conversions from [Jester et al.| 
l |2005[ and Lupton (2005, unpublished) which are available at 
https://www.sdss3.org/dr8/algorithms/sdssUBVRITransform.php. 
The zero points were then used for determining the magnitudes 
of the lensed quasars and lens galaxies. The results of the relative 
astrometry and photometry are given in Table The best fit 
parameters of the Sersic model assumed for the lens galaxies are 
reported in Table 
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Figure 1. Spectra of the 24 confirmed quasar pairs, including 13 confirmed gravitationally leased quasars (see Table In each panel, vertical dotted lines 
show the wavelengths of the quasar emission lines. Vertical short bars indicate the wavelengths of absorption features in the spectrum of the lens galaxy, which 
are shown only for lens systems for which lens galaxies are sufficiently bright to be detected in one of the quasar spectra. The redshifts shown in each panel 
con'espond to the quasar. The flux F\ is in units of ergcm“^s“^A“^. 


3.3 Notes on interesting individual systems 

Here, we discuss the properties of some of the interesting systems 
from the BQLS. 

3.3.1 SDSSJ0818+0601 

The spectrum shown in Figure[^indicates that the two stellar com¬ 
ponents separated by A6 ~ 1" have similar spectra (AV = 
226 km s“^), including broad absorption features in the Lya and 
CIV emission lines. However, no lens galaxy was found in our 
SOAR images (420 sec in i). If this is a gravitational lens, the 
non-detection of the lens suggests it is at a relatively high redshift 
(2i > !)• 

3.3.2 SDSSJ0821+0735 

While the spectra of the components (Figure [TJ indicate that the 
quasar redshifts are almost same {AV = 319 km s“^), their SEDs 
appear to be slightly different. The emission line features have dif¬ 
ferent profiles and the continuum shows differences in the slope. 
Some differences in the SEDs of lensed counterparts can be at¬ 
tributed to dust extinction in the lens galaxy (e.g., |Falco et al.| 


|1999| l or micro-lensing. However, the UH88 imaging observation 
(1440 sec in 1) also failed to detect any lens between the quasar 
components. 


3.3.3 SDSSJ0921+2854 

After we confirmed the lensing nature of this object, an archival 
search revealed that this lens was independently discovered by E. 
Ofek et al., and was observed by the Hubble Space Telescope Cy¬ 
cle 20 program GO-13001 although their results are unpublished. 
There are also X-ray and radio detections for this system, probably 
emanating from the lensed quasar. 


3.3.4 SDSS J0930+4614 

This is a quasar pair at the same redshift (AV =111 km s”'^) but 
with quite different SEDs (Eigure[T](. The fainter component has 
broad absorption lines associated with Lya, SilV, and CIV emis¬ 
sion lines, while the brighter component has no broad absorption 
features. Given the markedly different SEDs, this is likely to be a 
physical quasar pair rather than a lens. 
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Figure 2. GALFIT modelling results for the quasar lens sample using the reddest available band as labelled in each case. For each lens system, the columns 
from left to right are the data, the GALFIT model including the lens galaxy and the residuals with and without the lens galaxy, respectively. The flux scales 
in the data and model images match but they are different from the residual images. The flux scales of both residual images match. All images have the 
standard orientation with North up and East on the left. All images are 51 pixels on the side, which is ~n"for all systems except for JI442+4055 (~14") and 
J1458-0202 (~8"). 


3.3.5 SDSSJ1043+4320 

This is a quasar pair with similar redshifts (AV = 1778 km s“^). 
There are strong absorption lines in the CIV emission lines of both 
components, but the overall SED shapes appear to be different (see 
Figure[T](. It is important to conduct a deep imaging to search for a 
possible lens galaxy. 

3.3.6 SDSSJl 124+5710 

The SEDs of this quasar pair appear similar but have slightly differ¬ 
ent redshifts {AV = 496 km s“^, Figure[^. There is no indication 
of a lens galaxy in our UH88 imaging (1440 sec in /). Hence, we 
conclude that this is likely a binary quasar system. 

3.3.7 SDSSJl309+5617 

While the redshifts of the two components are almost identical 
(AV = 169 km s“^), the significantly different shapes of the Lyo? 


and CIV emission lines mean that this system is likely to be a bi¬ 
nary quasar system rather than a lens. 


3.3.8 SDSSJ1405+1350 

The SED of the fainter quasar has some similarity with the SED of 
the brighter quasar but the emission features have very low signal- 
to-noise ratios. The redshifts are not identical but imply that the 
quasars are probably physically associated (AV = 1400 km s“^). 
Additionally, the non-detection of a putative lens galaxy suggests 
this might be a binary quasar system. 
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Table 3. Relative astrometry and photometry of the 1 3 confirmed lens systems 


Object ID 

Component 

Ra 

Dec 

Red 

Blue 

and filters 


(J2000) 

(J2000) 

(7, i or R) 

y, g otR) 

SDSS JOl 14-t0722 

A 

01:14:38.421 

4-07:22:28.50 

20.88±0.07 

22.70±0.17 

ly 

B 

01:14:38.322 

4-07:22:29.35 

21.20±0.03 

22.98±0.09 


G 

01:14:38.386 

4-07:22:28.60 

19.30±0.03 

20.56±0.03 

SDSS J0256-t0153 

A 

02:56:40.771 

4-01:53:29.28 

21.20±0.05 

22.60±0.09 

ly 

B 

02:56:40.653 

4-01:53:30.06 

20.79±0.01 

22.33±0.04 


G 

02:56:40.740 

4-01:53:29.12 

19.89±0.04 

20.50±0.89 

SDSS J0737-t4825 

A 

07:37:08.677 

4-48:25:51.45 

18.28±0.00 

19.26±0.00 

ly 

B 

07:37:08.829 

4-48:25:51.26 

20.58±0.07 

22.03±0.14 


G 

07:37:08.782 

4-48:25:51.32 

19.87±0.10 

22.11±0.15 

SDSS J0821-t4542 

A 

08:21:58.706 

4-45:42:44.40 

20.76±0.14 

20.01±0.02 

Ry 

B 

08:21:58.582 

4-45:42:44.07 

20.86±0.02 

21.54±0.07 


G 

08:21:58.678 

4-45:42:44.38 

20.01±0.06 

- 

SDSS J0921-t2854 

A 

09:21:15.494 

4-28:54:44.57 

18.03±0.01 

20.36±0.01 

ly 

B 

09:21:15.355 

4-28:54:45.10 

18.93±0.01 

20.93±0.01 


G 

09:21:15.433 

4-28:54:44.35 

18.28±0.01 

20.28±0.08 

SDSS J1254-tl857 

A 

12:54:40.445 

4-18:57:12.33 

21.98±0.18 

_ 

R 

B 

12:54:40.295 

4-18:57:11.44 

21.74±0.07 

- 


G 

12:54:40.385 

4-18:57:11.89 

20.02±0.11 

- 

SDSS J1330-t3800 

A 

13:30:07.330 

4-38:00:43.82 

20.16±0.11 

19.52±0.02 

Ry 

B 

13:30:07.343 

4-38:00:42.39 

20.03±0.04 

19.69±0.02 


G 

13:30:07.375 

4-38:00:43.64 

19.74±0.15 

- 

SDSS J1412-t5204 

A 

14:12:10.158 

4-52:04:23.50 

18.91±0.01 

19.85±0.00 

ly 

B 

14:12:10.080 

4-52:04:22.16 

20.98±0.09 

22.00±0.02 


G 

14:12:10.114 

4-52:04:22.46 

20.16±0.06 

- 

SDSS J1442-t4055 

A 

14:42:54.795 

4-40:55:35.74 

17.64±0.00 

18.14±0.00 

i,g 

B 

14:42:54.614 

4-40:55:35.18 

18.28±0.01 

18.83±0.00 


G 

14:42:54.696 

4-40:55:35.38 

18.83±0.05 

20.07±0.10 

SDSS J1452-t4224 

A 

14:52:11.500 

4-42:24:29.02 

20.08±0.02 

21.05±0.02 

1,R 

B 

14:52:11.490 

4-42:24:30.61 

20.91±0.02 

21.71±0.02 


G 

14:52:11.490 

4-42:24:29.62 

18.15±0.02 

18.68±0.06 

SDSS J1458-0202 

A 

14:58:47.580 

-02:02:05.16 

21.50±0.02 

24.50±0.10 

i,g 

B 

14:58:47.670 

-02:02:03.49 

23.16±0.15 

22.76±0.02 


G 

14:58:47.650 

-02:02:04.19 

18.96±0.07 

21.51±0.08 

SDSS J1537-t3014 

A 

15:37:34.484 

4-30:14:53.70 

22.02±0.15 

21.95±0.10 

Ry 

B 

15:37:34.375 

4-30:14:55.49 

21.14±0.03 

21.97±0.18 


G 

15:37:34.442 

4-30:14:54.12 

18.95±0.29 

20.84±0.06 

SDSS J2146-0047 

A 

21:46:46.025 

-00:47:44.10 

19.74±0.02 

20.29±0.04 

ly 

B 

21:46:46.018 

-00:47:45.48 

20.54±0.05 

21.57±0.19 


G 

21:46:46.020 

-00:47:44.82 

21.37±0.32 

21.29±0.31 


3.3.9 SDSSJ1452+4224 

At Zs ~ 4.^ the quasar redshift of this lens is one of the highest 
redshifts known and is comparable to the lens SDSS J0946+1835 
( [McGreer et al.|20T0| ). The total lensing magnification is predicted 
to be ~ 8 by our mass models in Section]^ The lens galaxy is quite 
visible in the spectrum and the image, and has one of the lowest 
measured lens redshift {zi = 0.382) in our sample. Coinciden- 

^ Redshifts from Lya at such high redshifts are strongly affected by the 
Lya forest ^McGreer et al.|2010^ . Therefore, we give fewer significant dig¬ 
its on the redshift estimate. 


tally, the lens galaxy of SDSS J0946-I-1835 is at a similar redshift 
but more massive given the wider image separations of the lensed 
quasars. 


3.3.10 SDSS J2146-0047 

While the lensing nature of this object was confirmed by our 
imaging and spectroscopic follow-up observations, |Agnello et al.| 
( |2015b[ l recently reported an independent discovery of this lens us¬ 
ing Dark Energy Survey data. They tentatively assigned a lens red¬ 
shift of zi = 0.799 based on Mgll and Fell absorption features in 

























BOSS quasar lens survey 9 


the quasar spectra. The ahsorption feature is also seen in the Subaru 
spectrum shown in Figure[T] 

4 MASS MODELLING 

We created mass models for the 13 spectroscopically confirmed 
lens systems with the modelling software, GLAFIC ( |Oguri|2010| l. 
The redshifts of the quasars and the lens galaxies were taken from 
Table For those systems with no lens redshifts, we adopted a 
fiducial value of zi =0.5 (~ peak of the redshift distribution of 
typical lenses e.g., see Fig[^, as a lens redshift is required in the 
estimate of the velocity dispersion. 

We assumed that the mass distribution of the lens galaxies is 
isothermal. More specifically, we either used the singular isother¬ 
mal sphere (SIS) or ellipsoid (SIE) model for the lenses. Parameters 
such as the mass, ellipticity and the position angle (9e) of the lens 
model and the true (unlensed) position of the lensed quasar are al¬ 
lowed to be free. The positions and flux ratio^of the lensed quasar 
images served as constraints. We always begin with the SIS model 
and add further complexity such as adding more parameters only if 
the SIS results in a poor fit. The best fit parameters from the mass 
models are reported in Table[^ 

For SDSS J1330-I-3800, neither the SIS nor the SIE model 
could fit the constraints well. Hence, for this system, we used an 
SIS model with the lens galaxy position allowed to vary within the 
uncertainties from the GALFIT model. In a few other cases, we also 
used priors on the position angle of the lens based on the photo¬ 
metric models (see Table [^. In general, we could use the GALFIT 
errors on the positions and fluxes of the lensed images but had to 
relax thenj^ in a few cases in order to find a reasonable ■ The 
errors were relaxed up to ~ 0”15 and were always smaller than the 
pixel size. The best fit Einstein radii (or the velocity dispersions) of 
the BQLS lenses, which is a proxy for the mass of the lens galaxy, 
are similar to typical galaxy-scale lenses (e.g. |Bolton et al.|2006a^ . 

5 DISCUSSION 

5.1 Properties of the BQLS 

The basic properties of the 13 quasar lenses are summarized in Ta¬ 
ble We also compare the image separation, magnitude, and red¬ 
shift distributions of the newly discovered quasar lens systems with 
those of the SQLS sample in Figure]^ In the redshift comparison, 
we exclude systems for which the lens redshifts are not known. 
There are two SQLS lens samples ( |Inada et al.|2012| ): i) “all” - the 
entire lens sample which includes some heterogeneously selected 
lenses and ii) “stat” - this is a statistically well-defined sample with 
high completeness. The SQLS statistical lens sample is limited to 
the quasar redshift range of 0.6 < z < 2.2, i-band magnitude 
i < 19.1, image separations 1” < AO < 30”, and flux ratios for 
doubles larger than 10“° ®, in order to accurately characterize the 
selection function. 

Since the BOSS sample size is much smaller, we show nor¬ 
malised distributions to facilitate qualitative comparison of the 

^ The flux ratios are often known to be affected by dust, substructure, in¬ 
trinsic quasar variability and micro-lensing. Thus, a simple globally smooth 
models such as the ones tested here may fail to fit the flux ratios and cannot 
necessarily be used to rule out the lens hypothesis. 

This is reasonable because the errors from GALFIT are known to be an 
underestimate (e.g.,|Rusu et al.|2015| 


Table 5. Best fit lens mass model parameters. 


Object 

ID 

^Ein 

CO 

e 

Be 

(deg) 

SDSS JO 114+0722 

0.83 

0.42 

45 

SDSS 10256+0153 

1.03 

0.74 

42 

SDSS 10737+4825 

0.81 

0.25 

[-82] 

SDSS J0821+4542 

0.63 

0.51 

-21 

SDSS J0921+2854 

0.96 

0.53 

28 

SDSS 11254+1857 

1.15 

- 

- 

SDSS 11330+3800 

0.89 

[-0.553] 

[0.321]t 

SDSS J1412+5204 

0.84 

0.45 

52 

SDSS J1442+4055 

1.03 

- 

- 

SDSS 11452+4224 

0.87 

- 

- 

SDSS 11458-0202 

1.15 

0.27 

14 

SDSS J1537+3014 

1.14 

0.21 

[88] 

SDSS J2146-0047 

0.71 

- 

- 


The parameter ^Ein. e. Og denote the Einstein radius, ellipticity, and 
the position angle of the mass distribution, respectively. The numbers in 
bracket indicate which priors were used from the GALFIT model (see 
Table|^. t - This is an SIS model with priors on the lens position. Instead 
of the ellipticity and position angle, the last two columns list the best fit 
offset (in arcseconds) of the lens from one of the lensed quasars in x and y 
direction, respectively. We use the standard cartesian convention where x 
and y are positive to the right and up, respectively. 


Table 6. Summary of the confirmed quasar lenses. 


Object ID 


2; 

A6l(") 

i 

SDSS JOl 14+0722 

1.828 

0.408 

1.70 

19.1 

SDSS J0256+0153 

2.600 

0.603 

1.93 

20.0 

SDSS J0737+4825 

2.892 

- 

1.54 

18.5 

SDSS J0821+4542 

2.066 

0.349 

1.35 

18.9 

SDSS J0921+2854 

1.410 

0.445 

1.91 

18.9 

SDSS J1254+1857 

1.717 

0.555 

2.32 

19.5 

SDSS J1330+3800 

2.254 

- 

1.44 

20.1 

SDSS J1412+5204 

2.952 

- 

1.53 

19.0 

SDSS J1442+4055 

2.575 

- 

2.13 

17.9 

SDSS J1452+4224 

4.819 

0.382 

1.59 

19.1 

SDSS J1458-0202 

1.724 

- 

2.15 

19.1 

SDSS J1537+3014 

1.553 

0.490 

2.30 

19.6 

SDSS J2146-0047 

2.381 

0.799 

1.39 

20.2 


The image separations A9 are computed from Tablej^and the i-band PSF 
magnitudes are from the SDSS database. 

properties of these samples. The image separation distribution sug¬ 
gests that the BOSS sample finds a larger fraction of small im¬ 
age separation (AS < 2”) lenses compared to the SQLS sample. 
One likely reason for this difference is that we preferentially se¬ 
lected small-separation lens candidates for follow-up confirmation 
because larger separation quasar pairs are more likely to be physical 
pairs rather than real lenses (e.g.,|Kochanek et al.|1999]|Hennawi| 
|et al.|[2006a| ). Another possibility may be fiber collisions. Large- 
separation lenses are produced by clusters of galaxies typically 
at 2 ~ 0.5, and therefore there are usually many member galax¬ 
ies around the quasar images that are selected as Constant MASS 
(CMASS) galaxies for spectroscopy. This may reduce the chance of 
these quasar images being spectroscopically observed by BOSS. As 
mentioned earlier, this was indeed the case for SDSS J2222-F2745 
where none of the quasar images were observed by BOSS while 
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Table 7. Summary of confirmed quasar pairs. 


Object ID 

RAb 

Dccb 

ZA 

ZB 

A0(") 

iA 


A V(kms-Q 

SDSS J0213-0421 

02:13:22.865 

-04:21:34.341 

1.911 

0.992 

2.0 

20.8 

18.8 

> 2000 

SDSS J0818-t0601 

08:18:30.420 

+06:01:37.860 

2.359 

2.361 

1.1 

18.1 

20.0 

226 

SDSS J0821-t0735 

08:21:43.241 

+07:35:45.108 

2.383 

2.387 

1.3 

20.7 

21.6 

319 

SDSS J0928-t4332 

09:28:39.048 

+43:32:42.141 

3.694 

2.995 

1.6 

21.1 

20.3 

> 2000 

SDSS J0930-t46I4 

09:30:20.986 

+46:14:23.260 

2.393 

2.394 

1.5 

18.4 

19.9 

111 

SDSS J1043-t4320 

10:43:25.025 

+43:20:48.958 

2.245 

2.225 

1.7 

21.4 

20.7 

1778 

SDSS J1124-t5710 

11:24:55.440 

+57:10:58.120 

2.309 

2.315 

2.0 

18.6 

19.7 

496 

SDSS J1309-t5617 

13:09:27.338 

+56:17:41.161 

2.513 

2.515 

2.8 

20.5 

21.6 

169 

SDSS J1405-tl350 

14:05:56.865 

+13:50:39.827 

2.345 

2.361 

1.7 

21.7 

22.0 

1400 

SDSS J1548-t2830 

15:48:50.776 

+28:30:12.714 

3.208 

1.487 

1.7 

20.8 

20.9 

> 2000 

SDSS J2245-t2548 

22:45:55.801 

+25:48:33.548 

2.995 

2.178 

2.0 

20.4 

20.7 

> 2000 


See Table[^for the RA and Dec of the brighter quasar image (A). The positions and the i-band magnitudes in this table are from the SDSS database except 

for J0818+0601 where we used the astrometrically calibrated SOAR images. 





Quasar redshift (z^) 



Figure 3. Comparison of the BQLS with the SQLS “statistical” and “all” samples. Owing to the smaller sample size of BQLS, a qualitative comparison is 
more meaningful. Each histogram is thus normalised and integrates to unity. Lenses with no known lens redshifts are excluded from the panel showing the 
redshift distributions. BQLS lenses are fainter and have smaller image separations. The quasar redshift distribution and the peak of the lens redshift distribution 
of the BQLS sample are similar to the SQLS “all” sample. 


spectra were taken of several nearby member galaxies of this clus¬ 
ter lens. 

We find that 7 out of the 13 confirmed quasar lens systems are 
located at 2 : > 2.2. The fraction of high-redshift quasar lenses is not 
particularly high compared with the “all” SQLS sample (see Fig¬ 
ure]^. This result may be partly explained by the fact that some of 


the new lens systems were originally targeted as CMASS galaxies 
rather than quasars because of the dominant lens galaxy compo¬ 
nent. Indeed, for many of these cases, the lens galaxies are visible 
in the SDSS spectra which then provide the lens redshifts for these 
systems (see Table [TJ. Spectroscopic quasars selected in this man¬ 
ner are not necessarily high-redshift quasars. This also explains the 
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peak of the lens redshift distribution at zi ~ 0.5, which is close to 
the median redshift of the BOSS CMASS galaxy sample. 

Most of the BQLS lenses are fainter because BOSS target se¬ 
lection selects fainter quasars. The BOSS DR12 quasar catalogue 
contains 25 previously known lenses, many found by the SQLS. 
About two-thirds of this sample are bright lens systems (i < 18) 
with quasars at low redshifts (z < 2.2). We excluded this known 
sample before the visual selection of the BQLS lens candidates. 

A sub-sample of the BOSS quasars are selected homoge¬ 
neously and have a well-defined selection function which is called 
the CORE sample ( |Ross et al.|2012| l. Lenses present in the CORE 
sample could be used to define a statistical sample for the BQLS 
provided the visual selection process is calibrated. However, only 
two out of the 13 confirmed lenses are in the CORE sample (using 
the bossJargetl flag for QSOXOREJMAIN). It is unlikely that our 
visual inspection process is highly inefficient, so it is unlikely that 
we are missing a large number of true lenses. Therefore, even if we 
were to detect most of the true lenses, a statistically well-defined 
sample derived from this parent sample will still be probably too 
small for any statistical study. The “point-source” requirement for 
the CORE sample is probably the cause of this low yield. 


5.2 Quad fraction 

All of the 13 confirmed quasar lenses are two-image quasar lens 
systems. At galaxy scales, we expect about 15—20% of the lenses 
to be quadruply imaged for the BOSS quasar lens sample with 
f < 20 — 21 ( |Oguri & Marshall|2010^ which means we should ex¬ 
pect two or more quad lenses in our sample. There are several possi¬ 
ble explanations for the lack of quad lenses in the BOSS quasar lens 
sample. First of all, the BOSS quasar target selection selects only 
point sources to increase its efficiency ( |Ross et al.|2012| l, whereas 
small-separation quad lens systems are often classified as extended 
sources in the SDSS ( |Oguri et al.|2006| ). Indeed, all the four small- 
separation (AS < 2”) quad lenses discovered in the SQLS (see 
|Inada et al.|2012^ are found to be classified as extended sources in 
the SDSS dataset. In addition, the visual inspection method might 
be biased against selecting quad lenses because the quads appear 
to have odd morphology. And, spectroscopic follow-up observa¬ 
tions might be biased because single long-slit mode is used to tar¬ 
get the double quasar candidate images along with the candidate 
lens galaxy for confirmation. Therefore it is of great interest to ap¬ 
ply more sophisticated algorithms (e.g., |Chan et al.|2015[ |Ag nello| 
|et al.|2015a^ to the quasars in BOSS to search for quad lenses that 
are missing in the current lens sample. 


5.3 Quasar pairs 

Our follow-up spectroscopy identified an additional 11 quasar pairs 
with small angular separations of A9 < 3" (see Table for a 
summary). We expect most of these to be either projected pairs 
or binary quasars, based on their different SEDs and the lack of an 
obvious lens. Among the 11 quasar pairs, 7 pairs have the veloc¬ 
ity difference AV smaller than 2000 km s“^ and therefore satisfy 
the criterion of physical binaries used in |Hennawi et al.H2006a| l. 
The velocity difference may also be useful for distinguishing bi¬ 
nary quasars from lensed quasars. We find that the typical AV 
of the confirmed lenses measured from our follow-up spectra is 
< 40 km s“^ and ~ 100 km s“^ in a few cases, which can be re¬ 
garded as a typical error on our velocity difference measurements. 
From Table |7] we find that most of the quasar pairs have velocity 


differences significantly larger than those of the confirmed lenses, 
indicating that they are less likely to be lensed images of a single 
source but rather are physically associated distinct quasars. 

In general, there are still only a small number of quasar pairs 
known ([Hennawi et al.||2006b[ |2010[ [Myers et al.||2008] jKayoj 
|& Oguri|2012^ , and our sample adds significantly to the existing 
pair samples. For comparison, the SQLS only identified 8 binary 
quasars with A6 < 3" (see Kayo & Oguri|[20T2 for a compila¬ 
tion), which is comparable to the number of new binary quasars 
from the BQLS. The small-separation binary quasars have been 
used to study the very small scale clustering of quasars to discuss 
the possible role of mergers in enhancing quasar activities ([Hi 


et al.|2006b Hopkins et al.|2008}|Myers et al.|2008{|Shen et al. 


2010| Kayo & Og uripo" l. Projected quasar pairs are also useful 

for studying the distribution of absorbers around quasars ( |Hennawi| 
jet al.|2006a^ . 


6 SUMMARY AND CONCLUSION 

We present the initial results from BQLS, a systematic search 
for lensed quasars in the SDSS-III/BOSS data. We applied the 
same technique that we used for finding lensed quasars in the 
SDSS DR7 (SQLS). Here, we report the discovery of 13 confirmed 
quasar lenses. In addition, we present the discovery of 11 quasar 
pairs, some of which may consist of physical binaries. This sub¬ 
sample may still contain a few unrecognised lenses. The sample 
of new lenses includes one of the highest redshift lensed quasar 
(SDSS J1452-I-4224, z « 4.8) found to date. SDSS J1442-I-4055 
was also discovered independently by jSergeyev et al.] ( |2016[ . All 
the confirmed lenses from the BQLS have only two images. The 
lack of quad lenses is probably because BOSS only selects point¬ 
like quasars. We note that our follow-up observations are still in¬ 
complete with an additional ~ 50 good lens candidates that need 
verification. However, this sample is likely to have fewer real lenses 
since we conducted the follow-up observations of the most promis¬ 
ing candidates first. 

Compared to the SQLS, the selection function of the BQLS 
sample is complicated owing to the complex selection of the ini¬ 
tial BOSS quasar targets. In addition to the incomplete follow-up 
observations, a complex selection function prevents us from using 
the BQLS sample for statistical studies of either the lens popula¬ 
tion or cosmology (e.g., |Oguri et al.|2012[ l. The CORE sample from 
BOSS, which is a uniform and well-defined quasar sample could be 
used to produce a statistical BQLS sample, yielded only two lenses. 
This is probably because the quasars are required to be point-like 
in the CORE sample. 

A qualitative comparison of the initial BQLS sample with 
the SQLS shows that BQLS lenses are somewhat fainter and have 
smaller image separations, presumably due to the BOSS selection 
criteria. Furthermore, the BQLS and SQLS lensed quasars have 
similar redshift distributions in spite of the higher redshift selec¬ 
tion used for BOSS quasars. This is because many of the lenses in 
our sample were photometrically identified as galaxies before they 
were spectroscopically confirmed by BOSS to be quasars. In addi¬ 
tion, about half of the lensed quasars have redshifts lower than the 
quasar redshift range targeted by BOSS. These lenses would have 
been missed in the absence of BOSS CMASS galaxy spectroscopy. 

This study provides useful guidance for ongoing quasar lens 
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surveys in Hyper Suprime-Cartj^ and Dark Energy Survey jThe| 
Dark Energy Survey Collaboration]|2005^ in which even fainter 
quasar lenses will be discovered. 
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